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Definitic n of Contravariant Velocity Components 
Ching-mao Hung* 

NASA Ames Research Center, Moffett Field, CA 
Abstract 

This is an old issue in computational fluid dynamics (CFD). What is the so-called 
contravariant velocity or cor travariant velocity component? In the article, we review the 
basics of tensor analysis and give the contravariant velocity component a rigorous expla- 
nation. For a given coordinate system, there exist two uniquely determined sets of base 
vector systems - one is the c war i ant and another is the contravariant base vector system. 
The two base vector systems are reciprocal. The so-called contravariant velocity com- 
ponent is really the contravariant component of a velocity vector for a time-independent 
coordinate system, or the contravariant component of a relative velocity between fluid and 
coordinates, for a time-depeadent coordinate system. The contravariant velocity compo- 
nents are not physical quantities of the velocity vector. Their magnitudes, dimensions, and 
associated directions are controlled by their corresponding covariant base vectors. Several 
2-D linear examples and 2-E mass-conservation equation are used to illustrate the details 
of expressing a vector with respect to the covariant and contravariant base vector systems, 
respectively. 

I. Introduction 

Written in generalized curvilinear coordinates £ = y)-.V = v( x -> y) ^-D inviscid 

Navier-Stokes equations are, in a strong conservative form, 
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where, p, u, v, p, e are convex tional physical properties, J is transformation Jacobian and 

U = u + tyV, V = v *u + rj y v (2) 

Vinokur (Ref. 1) introduced t he strong conservative form of Navier-Stokes equations to CFD 
in 1974, and pointed out that (U, V) were the contravariant components of velocity vector. 
Compared to the Navier-Stokes equations written in Cartesian coordinates (x,y), except 
for the additional inclusion of the transformation Jacobian J, (U,V) are the replacements 
of velocities (u, v), and hence Steger (Ref. 2) in 1977 treated (U, V) as velocities and stated 

that (17, V) were 

“ the so-called contravariant velocities along the f and r/ coordinates.” 

Since then, the term “contravariant velocities” or “contravariant velocity components was 
widely used in computational fluid mechanics and heat transfer, in general without further 
statement or explanation, (for instances, Refs. 3 - 5). The term “contravariant velocity 
components” could be interpreted either as, like Vinokur (Refs. 1), “contravariant (velocity) 
components”, (i.e. contravariant components of velocity) which is strictly a mathematical 
definition or as, like Steger Ref. 2), “the components of contravariant velocity” in which 
the “contravariant velocity” can be loosely defined. In this article “contravariant velocity 
components” is treated as tie latter and is interchangeable with “contravariant velocities”. 

Following a similar definition of Ref. 2, while applied for body surface, Ref. 6 ex- 
plained that 

“ - - - the contravariant velocities are the decomposition of velocity vector into 
components along the £ coordinate line, U, and along the r? coordinate line, 7. 

There is some ambiguity in the above explanation. It had never explained the means and 
procedures of “decomposition” of velocity vector into components. It only described the 
directions the components associated but had never prescribed its corresponding base vec- 
tor system. One would intuitively interpret (U, V ) as the decomposed physical components 
of velocity vector along the directions of coordinate lines (£,??)• In contrast to Refs. 2 and 
6, Ref. 7 stated (for 3-D case) that 
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“The contravariant vei ocity components U, V, W are in directions normal to constant 
f , 77, C surfaces, respectively.” 

The above statement contradicts the statements from Refs. 2 and 6 about the directions 
of contravariant velocity components associated with a coordinate system. As we shall 
see, neither of the above two statements (from Refs. 6 and 7 ) is strictly correct. The 
contravariant velocity components are neither the physical components of a velocity vector 
nor in directions normal tc constant surfaces, respectively. The objective of this 

article is to clarify what are the contravariant velocity components and to further explain 
their physical implications. 

To avoid confusion, we will start with defining terms based on the algebra of vectors. 
The term component of a /ector will be defined. With the definition of covariant and 
contravariant base vector s /stems for given coordinates, the general mathematical term 
“contravariant component” will be discussed and its counter term “covariant component” 
will also be explained. For simplicity and clarity, we will use several 2-D linear cases as 
examples. Here in the 2 -D examples, the ( x\x 2 ) coordinates are interchangeable with 
(£,77) coordinates. 


II. Algebra of Vectors 

In an 71-dimensional vector space a set of n linearly independent vectors ei, e2, 

-, e n is called a base vector system. Any vector A in the space can be expressed as a unique 

combination of the base vectors A = Ciei + c 2 e 2 H 1- c„e n , or using the summation 

convention, A = c^e;. The vectors Ciei,c 2 e 2 , , c n e n are component vectors of vector 

A in the ei,e 2 , , and e n directions, respectively. The coefficients, c : ,c 2 , , and 

Cn , are called components of vector A in the e x ,e 2 , , and e n directions, respectively, 

with respect to the base (e. , e 2 , , ©n)- 

For example, in Cartesian coordinates, a vector A = 5i + 7j implies 

ex = i, e 2 = j, ci = 5, c 2 = 7. 

Expressing the same vector in terms of the base vector ((i - j)/2, j) so as A = 10(i - 
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j)/2 + 12j, then 


( 3 ) 


ei = (i - j)/2, e 2 = j, ci = 10, c 2 = 12, 
or in terms of the base vectors (2i, i + j) as A = — (2i) + 7 (i + j), then 

e x = 2i, e 2 = i + j, c x = -1, c 2 = 7. (4) 

As illustrated above, the ii.variant vector A can be expressed in different base vector 
systems which may or may not be orthogonal and normalized. While the “component 
vectors” are physical vectors the magnitude of the components of a vector not only depend 
on the directions but also depend on the magnitude (and the dimension, discussed later) 
of the base vector system. Obviously, without the prescription of “base vectors”, the 
“components” of a vector have no meaning. 

III. C mtravariant and Covariant Base Vectors 

Let x l be the coordinates of a point and r be a position vector. As shown in Fig. 1, 
a covariant base vector, g, = gr, is tangent to its corresponding coordinate line, x\ 
and a contravariant base vector, g* = W, is normal to the other coordinate lines x 3 
(or the surface of x l = constant), j ^ i ■ For any given coordinates x l , its contravariant 
and covariant base vectors, g* and g { are uniquely determined. They have the reciprocal 
relation, 


g* • g j = 




0 if i ^ j 

1 if i = j 


( 5 ) 


where 5j is the Kronecker delta. Based on the reciprocal relation, Eq. (5), we have 


ip-b = (no summation in i). 

& gtl cos (s l >&) 

If they are orthogonal, ther | g i |= ^ and g* and gl are in the same direction. Note 
that, the sketches in Fig. 1 (and Fig. 2) only give a general indication of the directions of 
g* and gi and do not imply any thing about their magnitudes. 


4 



Indeed, we have dr = g tdx* and the differential of arc length ds is determined from 
ds 2 = dr • dr; therefore 

ds 2 = g t • gj dx l dx j = gn dxddxd 

Here, the entity of gij thus defined is called the metric tensor, because of its being the 
metric of the space of the coordinates. Similarly, We can define g 3 g g - Ob\ iouslv, 
is conjugate or reciprocal tensor of c, f , i.e. 

9 ik 9kj = 

g l • g j = 9 lk Zk ■ g j = 9 th 9kj = ti]) 

The individual g ik and g kj are the components of two fundamental tensors. Let g = |ffy| 
and G = \g ij \ be determinants of gij and g ij , then ^/g = lfy/G where v fg is the Jaco- 
bian of coordinate transformation from old to new and VG is the Jacobian of coordinate 
transformation from new to old. Details of these relations can be found in regular Tensor 
Analysis” books, such as Rffs. 8 and 9. 

IV. Cont^ ariant and Covariant Components of a Vector 

A given vector A is an invariant, and can be expressed in contravariant and covariant 
base vector systems as 


A =- A : gi + A 2 g2 + A 3 g3 l-A n g n 

= Ajg 1 + A2g 2 + A 3 g 3 !-A n g n 

or, in the summation convention, 

A = A l g x = Ai gh ( 6 ) 

Here a super- script index is used for “contravariant” and a subscript index is used for 
“covariant”. According to he definitions in Sect. II, A 1 g 1 ,A 2 g 2 , ,A n g n are com- 
ponent vectors of vector A in the directions of gi,g 2 , ■ respectively. Simi- 
larly. Aig 1 . -4 2 g 2 , -.A n g n are component vectors of vector A in the directions of 
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gi g2 respectively. As defined in the previous Section, g; and g l are two 

uniquely defined base vector systems for a given coordinate system x l . (The terms covari- 
ant and contravariant will be explanded below.) Hence, we define 

A i as contravariant ccmponents of vector A in the directions of gi, with respect to 
the covariant base vector system (gi, g 2 i n)- 

Ai as covariant components of vector A in the directions of g l , with respect to the 
contravariant base vector system (g^g 2 , )• 

For simplicity, we call A 1 contravariant components of vector A and A l covariant com- 
ponents of vector A, with t xeir association with covariant and contravariant base vectors 
(respectively) implicitly implied. The reasons why A 1 are called contravariant and A { co- 
variant components of a vec tor are rooted in fundamentals of tensor theory (Refs. 8 and 
9). The essential property of a “vector” (or in general a tensor) is that its components 
follow the law of coordinate transformation. Setting the base vectors gi as the standard 
of comparison, if in a coordinate transformation, the components are transformed by the 
same matrix that described the transformation of the base vectors g t , the components are 
called “covariant” and if the components are transformed by the inverse transposed matrix 
of the the matrix that described the transformation of the base vectors gi, the components 
are called “contravariant”. The proof of that components Ai are “covariant” and A 1 are 
“contravariant” is beyond t ie scope of this article.) 

One can show that, from Eqs. (5) and (6), 

A 1 = A ■ g l and A{ = A • gi (7) 

Contravariant and covarianl components, A\ A l are related to each other 

A 1 = g %i Aj , and Ai = g t j A J 

To obtain “physical components” of vector A, one has to decompose the vector with respect 
to the normalized base vector, so as 

A =: 4*— = Ai— — ( summation on i). 

|gi| |g l l 
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( 8 ) 


Thus the physical componei ts A 1 along g* direction and A{ along g l direction are 

A* =| g t \A l = and A t =| g * M* = ( no sum on i). 

Here, the super-script index in component A 1 and subscript index in A t do not imply 
“contravariant” and “covarimt” as adopted in the article. They don t follow the law of 
coordinate transformation, hence they are not components of tensors. 

Eq. (7) is useful for evaluation but commonly leads to a mis- interpretation ( as in 
Ref. 7) that A 1 is the component of vector A in the direction of g* and Aj the components 
of vector A in the directior of g*. This mis-interpretation stems from the common use 
of a dot-product with a unit vector, so-called “projection”, to calculate the component 
of a vector in the direction of that unit vector. In projection, it implies that a vector is 
decomposed into two component vectors; one parallel to and another perpendicular to a. 
given direction. (This is different from the physical components mentioned above, as A 1 
calculated in Eq. (8).) In Ec. (7), neither the contravariant base vectors nor the covariant 
base vectors are normalized, so that Eq. (7) is not a process of projection. 


V. 2-D Examples and Continuity Equation 

Now in 2-D, let a position vector r(£, rj ) with coordinate system (£, 7?) , then the 
detailed formula, with relations to the Cartesian quantities, are 


A = r i + uj = A x gi -I- A 2 g 2 = Aig 1 + A2g 2 , 

g 1 = = £ x i + C y J, g 2 = Vr? = Tfci + Vyj 

dr dr 

gl = = *€ * + Vi J’ §2 = ^ = x v i + Vn J 

Vg = x^ur, - y$x v 

A 1 = (A • g *) = £ x u A A = (A • g ) = rjx u + Vy v 

A\ — (A • gi) = x^u + y$v, A 2 = (A ■ g 2 ) = x v u + y v v 

Here A 1 is the contravariant component of A in the direction of f and A 2 is the con- 
travariant component of A n the direction of rj, with respect to the base vector sj stem 



(§L \ A-i is the covariant, component of A in the direction of V£ and A 2 is the covariant 

component of A in the dim tion of Vr/, with respect to the base vector system (V<£, Vry). 

The relations of contravariant and covariant components to physical quantities can 
be illustrated in 2-D as Fig. 2. In Fig. 2, the vector A = OC is expressed by physical 
component vectors OB = A g r )gi = A 1 gi and OD = (A • g )g2 — A g 2 in the 
directions tangent to corresponding coordinate lines, (gi, g2), or component vectors OE = 
(A • g^g 1 = g 1 and OF = (A • g 2 )g 2 = A 2 g 2 in directions normal to the other 

coordinate lines, (g\ g 2 ), respectively. (Here, a segment with underline, such as OB, is a 
vector and without underline, such as OB, is a scalar.) In Fig. 2, OBCD is a parallelogram 
formed by BC || g 2 and CD || gi, and OECF is a parallelogram formed by EC || g 2 and 
CF || g 1 . Thus the physical components are, with respect to system (gi,g 2 ), 

OB =| gi | A 1 , OD =| g 2 | A 2 , 

and, with respect to system (g^g 2 )) 

OE =| g 1 \Ai, OF =| g 2 | A 2 . 


It is interesting to see the di Terence of physical component from the projection of vector in 
certain direction. In contract to OB and OD as the physical components, the projections 
of vector OC in directions cf covariant base vectors are 


OM = (A • gi)/|gi| = Ai/|gi|, ON = (A • g 2 )/|g2| = A 2 /|g 2 |. 

Similarly, the projections of vector OC in directions of contravariant base vectors 
are 


OG - (A • g 1 ) /jg 1 1 = A7|g x |, OH = (A ■ g 2 )/|g 2 | = A 2 /|g 2 |. (9) 

We would like to illustrate the above discussions with several examples, by using the 
same vector as in Sect. II, A = 5i + 7j, expressed in various contravariant and covariant 

base vectors. 
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As a first example, consider the coordinate system £ — 2x , r/ — y which is still 
orthogonal but stretched, thus 

r(t,v) = ^/ 2 i + ??j 

g 1 = = 2i, g = Vx'i + byj = j 

gi = i + Vi j = l/2i, g 2 = x r ) i + y v j = j 
Therefore, the contravarianl components for the vector A are 

A 1 = A • g 1 = 10, A 2 = A-g 2 = 7 

and the covariant components are 

A L = A-gi = 5/2, A 2 = A • g2 = 7. 

Fig. 3 shows the sketch of tl is example, with magnitude and direction. Note that the new 
coordinates (£, rj) are same as the (x, y ) Cartesian coordinates, except £ is scaled as 2x. 
While g 2 and g 2 are identical with in both magnitude and direction, gi and g 1 are in the 
same direction as x-coordinite, i, but different in magnitude, and the magnitudes of A 1 
and A i are scaled accordingly. The shaded area in Fig. 3 (also in Figs. 4a and 5) indicates 
as the value of y/g, Jacobiar of transformation from old to new coordinates. 

As a second example, consider the coordinate lines (f , rj) defined by £ = 2x, r) = x+y, 

so as 


r (£> v) = £/2i + (^ ~ C/2)j 
g 1 = 2 i, g 2 = (i+j) 
gi = (i - j)/2, g2 = j 

Therefore, the contravariant components for the vector A are 

A 1 = 10, A 2 = 12 
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and the covariant components are 


A\ = — 1, A 2 7. 

These are identical as Eqs. 3) and (4) in Sect. 1, written as 

A = A l gi + A 2 g2 = 10(i — j)/2 + 12j 
A = J 4 1 g 1 + J 4 2 g 2 = -(2i) + 7(i + j) 

Fig. 4 shows detailed sketc res of the base vectors, with magnitude and direction, and 
relations with ( x , y). Fig. 5 shows the vector A expressed in its contravariant and covariant 
base vector systems. Since tie magnitudes of base vectors are not normalized, the marked 
number are scaled to each brse vectors accordingly. 

It is interesting to note that if vector A is decomposed into two component vectors 
so that one is parallel to g 1 and another is normal to g 1 . then A = 5/ \/2 + 7/2 ^ Jijp ■ 

These are the projection of ector A on to unit vectors (i — j)/ \/2 and (i -f j)/ \/2 , and are 
quite different from the expression in contravariant or covariant base vector system. 

As a third example, lei a new coordinate system be £ = (x — y)/ 2, rj = y, and hence 
r(£, rj) = (2/ + rj)i -|- r/j. Therefore, the contravariant and covariant base vectors are 

g 1 = (i-j)A g 2 = J 

gi = 2i, g 2 = (i + j) 

Hence, again as Eqs. (3) and (4), 

A = Aig 1 + A 2 g 2 = 10(i - j)/2 + 12j 
A = A 1 g 1 + A 2 g 2 = — (2i) + 7(i+j) 

Fig. 6 shows the contravariant and covariant base vectors for this new coordinate system. 
Note that the coordinate system in Fig. 6 is the conjugate or reciprocal of the the coor- 
dinate system in Fig. 4. The contravariant base vectors in the 3rd example are covariant 
base vectors in the 2nd example and the covariant base vectors in the 3rd example are 
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contravariant base vectors in the 2nd example. Clearly, the vector A is an invariant, and 
the specification of coordinate system is essential for vector decomposition. 

In fluid dynamics, the continuity equation can be written in integral and differential 
forms as 

^ f pdv+ [ pW • ndS = 0, 
dt J v Js 

J !/-(«)* + JjPW - ( 5 ) V + JCw ■ (g)h 2 = 0. < 10 ) 

Here, p is density, W = ui + nj the velocity vector and the contravariant velocity com- 
ponents, U = W • g 1 and / = W • g 2 are contravariant components of velocity vector 
W, in the directions of covariant base vectors gi and g 2 , respectively. Here y/g is the 
same as J~ l in Eq. (1), tht Jacobian of transformation from old to new coordinates (or 
g ~ l = y/G = J). In a finite volume formation, y/g is the volume of the grid cell and 
( 0)5 g 1 and (g)^ g 2 are surface normal in the directions of g 1 and g 2 , respectively. The 
concept of mass conservation can be interpreted accordingly. While ( U , \ ) are computed 
by the dot-product of velocity with g 1 and g 2 , respectively, it should now be clear that 
they are neither physical components of velocity vector nor in directions normal to the 
corresponding surfaces. 

It is important to note that the contravariant velocity components, (U,V), do not 
have the physical dimension of W (velocity) itself. In fact, in addition to their magnitudes, 
the dimensions of contravariant and covariant components are controlled by the dimensions 
of the covariant and contravariant base vectors, respectively. In Eq. (5) the contravariant 
and covariant base vectors are reciprocal not only in magnitude but also in dimension. This 
can be clearly illustrated by writing the continuity equation in cylindrical coordinates as 

£- r= yfix* + y 2 ), rj = 9 = tan~ 1 y/x 

gi = g 1 = cos9 i + sinQ j, yfg = r, 
g 2 = — r (sin9 i — cosd j ), g 2 = — 1/r ( sin9 i — cosO j ) 
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U = u cosd 4- v sinO , — \/r{y cosO — u sinO) 

d c) d 

— rp + —ro(u cosO + v sinO ) + — p(n cosO — u sinO ) = 0. 
at or oO 

Here, gi and g 1 are identical. Even they are same in the direction, (orthogonal to gi 
), g 2 and g 2 are different in magnitude and dimension. Obviously, in the above, the 
dimension of contravariant velocity component U is different from that of contravariant. 
velocity component V. 

VI. Time-dependent Coordinates and Different Variations 

Up to this point, we have only dealt with the cases where the coordinate system 
independent of time. Actually, Steger (Ref. 1) and Refs. 5 and 6 did include a general 
transformation with time at, in 2-D, £ = £(x, y , t), y = r?(: r, y, f), r = t and defined 
the contravariant velocity components as 


U — C "H ff 9i /V — 9.x ) d” Vt )i 

V = Tjt + rj x u + TjyV = T] x (u - x T ) + r/y(v - y T )■ 

The above relations can be rewritten as 

U = W-g 1 - w 7 -g 1 = !W- w') -g 1 , V = W • g 2 — w' • g 2 = (W — w') • g 2 (11) 

where w' = ^ is the coordirate velocity vector. The additional terms are the contravariant 
components of the coordinate velocity vector. Indeed, — w' is the relative velocity 
vector. Therefore, the so-called contravariant velocities are the contravariant components 
of relative velocities, between the fluid and the coordinates (or locally, the grid lines). 
Obviously, all the previous discussions still hold for the time-dependent coordinate system. 

Finally, we would like to point out that not all the so-called contravariant velocities 
are defined in the same form as discussed above. For instance, in Refs. 10 and 11, they 
were called as “contravariant velocities written without normalization”. Refs. 10 and 11 
have never clarified what would be the “normalized” contravariant velocities. In Ref. 12 
the contravariant velocities were defined as 

rr 1 cr - 

u = V = W-J~ (12) 

ig 1 ! |g~l 
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In Ref. 13 the “physical con ;ravariant velocity components" (£7, R) were defined as 

U = W-| gl |g\ R = W ■ jg2 |g 2 (13) 

In Refs. 14 and 15 by including the Jacobian of transformation, the contravariant velocities 
were defined as 

U = W-(g)*g\ R = W • {g)?g 2 (14) 

The three above definitions are quite different from each other. In Eq. (12), (U, R) are 
the projections of velocity vector on the directions of g 1 and g 2 , respectively, and are the 
same as OG and OH in Fig. 2 given by Eq. (9). They are not “components” of the velocity 
vector W as pointed out previously. In Eq. (13), (17, R) are the physical components 
of velocity vector along the f and r) directions, given by Eq. (8), and as OB and OD in 
Fig. 2, respectively. The (17 R) defined in Eq. (14) are not velocity components either. 
In a finite-volume formulation, Eq. (14) represent the velocity dot-product with surface 
normals (in the direction of the contravariant base vectors), i.e. the convective flux rate 
across the corresponding coitrol surfaces. The dimensions of U and R are still controlled 
by the base vectors g 1 and g 2 respectively. This can be seen in the continuity equation, 
Eq.(10). Nevertheless, the ibove three variations of (U, V) are no longer components 
of a tensor, and so are neither contravariant nor covariant components of a vector. We 
would rather hold to the original meaning of “contravariant” and maintain the definition 
of contravariant velocities as Eq. (2) or Eq. (11), without any modification. 

VII. Concluding Remarks 

In this paper we have reviewed the basics of tensor analysis in an attempt to clarify 
some misconception regard c mtravariant and covariant components of a vector as used in 
fluid dynamics. We have inc.icated that, for any given coordinate system, there are two 
uniquely determined reciprocal covariant and contravariant base vector systems. The con- 
travariant components are components of a given vector expressed as a unique combination 
of the covariant base vector sj stem and vis versa, the covariant components are components 
of a vector expressed with the contravariant base vector system. A vector is an invari- 
ant. It is the components which are contravariant or covariant, associated with the given 
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coordinate system. And the so-called “contravariant velocity” or “contravariant velocity 
component” is really the contravariant component of the velocity vector. As described in 
Section of Algebra of Vectors, expressing a vector with a combination of base vector is a 
decomposition process for a specific base vector system. Hence, the contravariant velocity 
components are decomposer components of velocity vector along the directions of coordi- 
nate lines, with respect to the covariant base vector system. However, the contravariant 
(and covariant) components are not physical quantities. Their magnitudes and dimensions 
are controlled by their corresponding covariant (and contravariant) base vectors. 

For a time-dependent coordinate system, the so-called contravariant velocity com- 
ponents are the contravariant components of the velocity vector minus the contravariant 
components of the coordina te velocity vector. In other words, the contravariant velocities 
can be interpreted as the contravariant components of relative velocities, between the fluid 
and the coordinates (or loca ly, the grid lines) . Obviously, all the discussions hold for both 
time-dependent and indepei dent coordinate systems. 
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Fig. 1 Sketch of Cont^avariant and Covariant Base Vectors, g 1 and g 
( indication of direction only). 



Fig. 2 Relations of Physical Components and Projections to 


Contravariant and Covariant Components 



A = OC = OB + OD = OE + OF 
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